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CN . Abstract 

(N 

^D ' The zero neutrino mode of the double beta decay in heavy deformed nuclei is investigated 

in the framework of the pseudo SU(3) model, which has provided an accurate description of 

■^ i collective nuclear structure and predicted half-lives for the two neutrino mode in good agreement 

0^ ' with experiments. In the case of 238[/ ^j^g calculated zero neutrino half- life is at least three orders 

of magnitude greater than the two neutrino one, giving strong support of the identification of 

the radiochemically determined half-life as being the two neutrino double beta decay. For 

CJ I ^^^Nd the zero neutrino matrix elements are of the order of magnitude of, but lesser than, those 

H ■ evaluated using the QRPA. This result confirms that different nuclear models produce similar 

zero neutrino matrix elements, contrary to the two neutrino case. Using these pseudo SU(3) 

results and the upper limit for the neutrino mass we estimate the /3/3o!^ half-lives for six nuclei. 

/\f • An upper limit for majoron coupling constant is extracted from the experimental data. 



S 



1 Introduction 



Many grand unified theories |]|, ^, ^ regard neutrinos as Major ana particles, identical to its 
own antiparticle, predicting light neutrino masses and the existence of right-handed currents. 
In order to get information about this physics beyond the standard model many extensive 
experimental investigations have been performed to detect the neutrinoless mode of the double 
beta decay (/3/3oiy)- 

Double beta decay can be classified into various modes according to the light particles besides 
the electrons associated with the decay. The two neutrino mode {(3(32v)-, in which two electrons 
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and two neutrinos are emitted, takes place independently of the neutrino properties, and it is 
the only one which has been measured up to now @, ^, y, 0, g, y, |l^ . 

The nonconserving lepton number zero neutrino mode (/3/3ojy), in which only two electrons 
are emitted, can occur if, and only if, the neutrino is a massive Majorana particle W^, i-e. if 
u and P are just the left and right handed helicity states of a single particle state i/^^For this 
reason this mode, unobserved up to now, concentrates the interest concerned with new physics, 
beyond the standard model. 

Another proposed double beta decay mode is the decay in two electrons and a massless 
neutral pseudoscalar X) known as the majoron. If this PPoux decay competes with the PP2u 
mode, the electron energy spectra of the two emitted electrons must be compared in detail in 
order to distinguish both decay modes. 

In previous papers [^] we used the pseudo SU(3) shell model to evaluate the two neutrino 
double beta half lifes of eleven heavy deformed potential double beta emitters. We found good 
agreement with the available experimental information. The radiochemically measured ^38^ 
decay jl^ raised some expectatives, given the experimental array cannot discriminate between 
the different double beta decay modes, and some theoretical estimations [^] predicted similar 
decay ratios for both the Ov and the 2^ modes. Our calculations for the two neutrino mode 
were consistent with the experimental result. In the present work we estimate the zero neutrino 
matrix elements for six heavy deformed double beta emitters, including 238^/ g^^^^ ^^^Nd. In the 
case of 2^^f7 we found the zero neutrino half- life at least three orders of magnitude greater than 
the two neutrino one, giving strong support of the identification of the observed half-life as being 
the two neutrino double beta decay. Experiments in ^^^Nd have set the lowest limits to the 
majoron coupling constant. For ^^^Nd the f3f3oi, matrix elements are of the order of magnitude 
of, but lesser than, those evaluated using the QRPA. This result confirms the experience that 
different nuclear models produce similar zero neutrino matrix elements, contrary to the two 
neutrino case [^, |3|, |l^, 14, |l^. Using these pseudo SU(3) results and the upper limit for the 



neutrino mass we estimate the PPou half-lives for six nuclei. An upper limit for the majoron 
coupling constant is also extracted from the experimental data. 

In section 2 we review the double beta decay formalism for the zero neutrino and majoron 
emitting modes. A summary of the pseudo SU{3) model is given in Section 3 and in Section 

4 the calculations of the nuclear matrix elements Mqj^ within this model is presented. Section 

5 contains the results and discussions and Section 6 the conclusions. In Appendix A some 
details of the evaluation of the radial integrals are given, the connection with the 2u formalism 
is presented in Appendix B. The inclusion of short range correlations and the effects related 
with the finite nucleon size are discussed in Appendix C. 



2 Double beta decay 



In the case the Ov decay exists, the virtual neutrino must be emitted in one vertex, and absorbed 
in the other. Since in the standard theory the emitted particle is a right-handed antineutrino 
and the absorbed one a left-handed neutrino the process requires that a) the exchanged neutrino 
is a Majorana particle and b) both neutrinos have a common helicity component. The helicity 
matching can be satisfied in two ways: a) the neutrinos have a nonvanishing mass and therefore 
a "wrong" helicity component proportional to rriy/Ej^. The decay rate will be proportional to 
< rriy >^. Or b) the helicity restriction could be satisfied if there is a right handed current 
interaction. In this case a nonvanishing mass allowing mixing of neutrino types is also required 



©H]- 



For massive Majorana neutrinos one can perforin the integration over the four-momentum 
of the exchanged particle and obtain a "neutrino potential" which for a light neutrino (mjy < 
WAIeV) has the form 

H{r,E) = — / dq -^- (1) 

Trr Jo q + E 

where E is the average excitation energy of the intermediate odd-odd nucleus and the nuclear 
radius R has been added to make the neutrino potential dimensionless. In the zero neutrino 
case this closure approximation is well justified |]l^]. The final formula, restricted to the term 
proportional to the neutrino mass, is g, |l^] 

(ro?r^ = (^^)'Go.MoV (2) 

where Gq^ is the phase space integral associated with the emission of the two electrons and the 
nuclear matrix elements Mq^ are discussed below. 

The double beta decay with majoron emission has been investigated in various models||3|, pi] ]. 
The original Gelmini-Roncadelli one|l^ is incompatible with the LEP measurement of the 
invisible width of the Z boson. However, various suitable modifications have been considered [16, 



19|. The decay rate for the majoron emitting mode of the double beta decay resembles the Qv 



mode: 
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«;)"' =< 9.,x > Go.,xMo^.- (3) 

where < g^^-^ > is the effective majoron neutrino coupling constant, Gq^,^ is the phase space in- 
tegral which describes two electrons plus the massless majoron, and the nuclear matrix elements 
are the same in both cases. 



These nuclear matrix elements are |11| 



Mo. ^ |Afo^/ - 4<J (4) 

9a 

with 

Mo^=<0+||0-||0+> (5) 

where the kets |0^) and |0t) denote the corresponding initial and final nuclear states, the 
quantities gy and gA are the dimensionless coupling constants of the vector and axial vector 
nuclear currents, and 

O^"" = E Ogl = E ffmt;n ■ <^nt-H{\fm " fn\,E) 

m,n m,n /p\ 

o^^E o^L = E t;;,t-H{\rm - fn\,E). ^^^ 

m,n m,n 

being a the Pauli matrices related with the spin operator and t~ the isospin lowering operator, 
which satisfies t~\n >= \p >. The superindex GT denotes the Gamow- Teller spin-isospin 
transfer channel, while the F indicates the Fermi isospin one. 

It must be noted that only the difference between Fermi and Gamow- Teller matrix elements 
appears in the combination Mq^, which can be evaluated from the new operator 



O = Y.t-,t-H{\r^ - rn\,E){am ■ < - {^—f) (7) 

m,n "^ 

Although (— )^ = 1.5 for free nucleons [O, 14], it is common to assume (— )^ = 1-0 in 
order to improve the agreement with the observed Gamow- Teller transition strength, which is 
quenched by roughly a factor one half in comparison with the theoretical sum rule [^ . In this 
way, channels which are out of the most usual models, as the delta-isobar-hole excitations, are 
included in an effective form |21]. Under this assumption 

O = E tmtnH{\rm -r„|,S)(<Tm ' o'n " 1) 

■>n,n ,g\ 

= -4 E t^t-H{\rr^ - fr,\,E)Pmn{S = 0) ^ ' 

m,n 

where PmniS = 0) is the projector of a two-particle state coupled to spin S=0. This result will 
be useful in order to simplify the expressions for the nuclear matrix element Mq^. 

The finite nucleon size and the two-nucleon short range correlations in the nucleon-nucleon 
interaction are introduced in these matrix elements. Their explicit expression are given in 
Appendix C. 

As mentioned above, we will use the pseudo SU(3) scheme to describe the ground state of 
the six heavy deformed nuclei considered in this contribution. A detailed description of this 



formalism can be found elsewhere [12, ^, 23]. In the next section we give a brief summary of 
the model. 



3 The Pseudo SU(3) Model 



A tractable shell model theory for deformed nuclei requires a severe truncation of the spher- 
ical model basis. A successful truncation can be achieved if the basis selection is made rela- 
tive to those parts of the interaction that dominate the low-energy structure: the pairing and 
quadrupole-quadrupole ones. In the case of rotational nuclei, it has been shown that the long 
range part , i.e. the quadrupole-quadrupole interaction, is the most important one [p^]. 

For heavy nuclei, the major shells are built by the orbitals j = {1/2, 3/2, . . . ry — 1/2}, which 
are called of normal parity plus a single particle state j = rj + 3/2 of a different parity. Then 
we are restricting the Hilbert space to the valence neutron and proton shells constituted by the 
corresponding normal parity orbitals and the abnormal or unique parity levels. 



In the pseudo SU{3) shell model coupling scheme ]22] the normal parity orbitals are identified 
with orbitals of a harmonic oscillator of one quanta less t/ = r] — 1. This set of orbitals, 
with j = j = I + s, pseudo spin s = 1/2, and pseudo angular momentum taking the values 
1 = 7], ry — 2, ... 1 or define the so called pseudo space and only recently it was found an 
analytic expression for the transformation that take us from the normal parity orbitals to the 
pseudo space |2J]. Applying this transformation to the spherical Nilsson hamiltonian it can be 
shown explicitly that the strength of the pseudo spin orbit interaction is almost zero and the 
orbitals j = / =b 1/2 are nearly degenerate doublets. For the configurations of identical particles 
occupying a single j orbital of abnormal parity, a convenient characterization of states is made 
by means of the seniority coupling scheme. 

The many particle states of n^ nucleons in a given shell 77^, a = i/ or vr , can be defined 
by the totally antisymmetric irreducible representations {l"" } and {1""} of unitary groups 



of dimensions corresponding to the sizes of the normal Q^ = {fja + !)(??« + 2) and unique 
O^ = 2ria; + 4 parity spaces, respectively; with the constraint n^ = n^ + n^ ■ 

A complete classification of the states can be defined by the following chains of groups: 

{1"^} {fa} {faha{Xa,f^a) ~Sa K^ la J^ 

Uin^) D U{n^/2) X U{2) D SU{3) X SU{2) D 50(3) x SU{2) D SUj{2) 

(9) 



{1"^} ^a 13a J^ 

u{ni) D Sp{ni) D su{2) ' 

where above each group the quantum numbers that characterize its irreducible representations 
(irreps) are given and the 7^, /?„ and K^ are multiplicity labels of the indicated reductions. 
For the normal parity spaces the pseudo LS coupling scheme is used and the following relations 
between its quantum numbers are satisfied 

< = fa + fl S^ = l/2{fl - fl\ {U} = {2</2-^^ l^s^} , (10) 

while for the abnormal parity spaces the seniority configurations Va are appropriate. 

Proton and neutron states are coupled in angular momentum in both the normal and unique 
parity sectors, generating states with angular momentum J and J respectively. The wave 
function of the many-particle state with angular momentum J and projection M is expressed 
as a direct product of the normal and unique parity ones, as: 

\JM>= J2 ilJ'' > ^IJ"" >]m (11) 

We will be interested in the ground states of some deformed nuclei, which have J = M = 0. 

For even-even heavy nuclei, it has been shown that if the residual neutron-proton interaction 
is of the quadrupole type, independently of the interaction in the proton and neutron spaces, 
for yrast states below the back bending region the most important normal parity configurations 
are those with highest spatial symmetry {fa} = {2""'^}. This implies that S-,^ = Sy = 0, that 
is, only pseudo spin zero configurations are taken into account. 

Additionally in the abnormal parity space only seniority zero configurations, w,r = v^ = 0, 
are taken into account. This simplification implies that (^^^ = (3^ = 1 and J^ = J^ = 0. This is 
a very strong assumption, which in future works is expected to be improved, but is quite useful 
in order to simplify the calculations. 

As an example, the ground-state band of 238[/ ig characterized by the eigenstates 

I {l6},{23},(18,0)^; {li2},{26},(36,0),; 1(54, 0)E:JM >^ 

I (ii3/2)^, J^ = M^ = 0; (ii5/2)^, J^ = M^ = 0>A , (12) 

where the first term indicates the normal part and the second the abnormal one. The ^38^ 
ground state has K = 1,J = M = 0. For a detailed description of the way this state was 
constructed the reader is referred to [^, |2^ . 

4 Calculation of the matrix elements Mqjj 

We want to evaluate the nuclear matrix element (4) for several heavy deformed nuclei, which are 
potential Pf3 emitters. According to the above discussion the initial [0^) and final lOt) ground 
states are of the type indicated in the Eq.(12). 



We expand the two body transition operator in order to obtain its normal and abnormal 
parity components 

O" = - En,j, E,3,i4 ^J,M < {jiJ2)JM\Of,\{hJ4)JM > [ajM'^naj,a,y^ 
= - Ej.j^, E>j,, Ejm < {JJ.')JM\Of,\Uuju')JM > [alal,Y"'[a,^a,J-J'^ 

where in the second line only the nonvanishing contributions were retained. Transforming this 
operator to the pseudo SU{2>) space, we have formally the expression 

O" = O^.^. + O^^^^ + OXj,^ + OXa^ (14) 

where the subindices A^A^, A^^, . . . are indicating the normal or abnormal spaces of the fermion 
creation and annihilation operators, respectively. By using the Nilsson scheme in order to obtain 
the occupation number we are grouping nucleons in pairs. For this reason only these four cases 
can give contribution different from zero. 

In the present work we will restrict our analysis to six potential double beta emitters which, 
within the approximations of the pseudo SU(3) scheme, are also decaying via the 2v mode. 
They include the observed "^^^Nd^ ^^^Sm and ^^^f/ —f '^^^Pu decays. In this case two neutrons 
belonging to a normal parity orbital decay in two protons belonging to an abnormal parity one. 
Under the seniority zero assumption for nucleons in abnormal parity orbitals, only proton pairs 
coupled to J=0 are allowed to exist in the ground state, restricting the above summation to 
J=0, M=0. It follows that vr = vr' and u = u' . The matrix elements becomes 

Mo^= - E <(i.i.)J = M = 0|O"|(i,>)J = M = 0> 

1 {2</2}(At,^^);{2"^/2},(AJ„/xt);(A/,/x/)100 > (15) 



In these cases, it is also simple to show that [12 



< U^r^^^ J^ = M^ = |[at at y=^^=^\ (j^r^, J^ = M^ = > 






^^ / («^+2)(2v+l-n^) 



(16) 



The operators in the normal parity sector must be transformed to the pseudo-space. These 
nucleons are coupled to pseudospin 5 = 0, which together with the J=0 condition implies L = 0. 
Then 111, H 



Expanding the normal states in their proton and neutron subspaces, and the above operator 



in their SU(3) components, together with Eq. (16), we arrive to the final expression 






E '"e < (A/,/x/)ir^L^,(A^,/x/)E::,L,||(A/,/i/)10 > 

< {Xi,f,i)K^L^,{Xi,^il)K,L^\\{X\fi')W > 

{\o,^J.o)po 

< (0,r?^)l[^,(0,r/^)l/;||(Ao,/io)10 > 

(18) 

where the < {Xi, fii)KiLi, {X2, fJ-2)K2L2\\{X, fi)KL >p are SU(3) Clebsch-Gordan coefficients, 
and < (Ai,/ii)|||0|[|(A2,/U2) > are triple reduced matrix elements [|^]. We have implicitly 
defined Mg^(jj^) as the contribution of each normal parity neutron state ji, to the nuclear 
matrix element in the transition (j^)^ — > {j^Y ■ 

The two-body matrix element can be expanded in its L, S components 

< {inJn)J = M = OiO"|(i .j.)J = M = 0> 

= E5 V(2i. + l)(2i. + 1)(25 + l)T^(ij^5/^; /.i)l^(ij.5/,; Z^) (19) 

< {IJ^)S\\H{t, E)\\{IJ,)S X (ii)S||r ■ r(a)||(ii)5 > 

where the W (...;...) are Racah coefficients and 

r • T{GT) = ai-a2 T- T{F) = 1 a = GT or F. (20) 

In order to evaluate the spatial matrix elements, we use the Bessel- Fourier expansion of the 



potential |25|, which gives 



<{hl2)LM\V{r)\{kk)LM> 

= E/(-l)'^+^*+^(^l||C/||/3)(/2||Q||M^(W2y4;W)(2^ + l)^'ai^2,/3/4) ^ ^ 

with (/1IIQII/3) the reduced matrix elements of the unnormalized spherical harmonics Cimi^) = 

v 2lTT^'™(^) ^^'^ R (hh^hh) the radial integrals which are explained in Appendix A. 

2 
Although we will present results for MqI and Mq^ , we have shown that using -j^ = 1.0 

only the 5 = channel survive, i.e. only states with L=0 are of interest and, after some 
manipulations, we obtain the simplified expression 

< iJnjn)J = M = 0|0|(>>)J = M = > 



(2i^+l)(2>+l) 



= -2 f(2t+i)(2i+i) < (^-^-)^ = mir,E)\iUu)L = > (22) 

= -2v/(2i, + l)(2j. + 1) Ei-'^yiluOLO\lO)^R\U^, IJ,). 

I 

The radial integrals R {It^Ij^jIi^I,^) are evaluated in the way described in Appendix A, and 
they include finite nucleon size and short range correlations, as explained in Appendix C. 



5 Results and discussion 



We evaluated the matrix elements of Mq^ between the ground states of six different initial 
and final nuclei for the /3/3 transitions. There exists a relevant dependence of this transition 
amplitudes with the partition {n^,n^ ,n^,n^} of the valence nucleons of each pair of nuclei. 
We used the standard procedure to fix them, taking deformation from experiment, and filling 



each Nilsson level with a pair of particles in order of increasing energies |12, ^. As mentioned 
in the Introduction, we restricted the present analysis to nuclei in which the unique and normal 
proton and neutron valence occupations are related by 

The six nuclei fulfilling these conditions are listed in Table 1, and are the same for which finite 
half-lifes were predicted in the two neutrino case [|l^]. Their matrix elements Mq^, are shown 
in the second column. The value for the integrated kinematical factor, Gqi, in [years" ^] was 



obtained following the procedure indicated by Doi et al |26] with qa/ Qv = 1-0 and it is given in 
the third column. In the fourth column of Table 1, the energy value E is indicated and in the 
last two columns the theoretical predictions ( for < rUe >= leV) and experimental lower limits 
of the /3;9ojy-half lives are given. It has been demonstrated that these Mq^, matrix elements are 
nearly insensitive to the exact values of the average energies if they are in the range of some 
MeV 0. We choose the accepted parametrization E = 1.12^4^/^ [||]. 

As mentioned above, for ^38^ ^^^ predicted Ov half live Tq[, = 1.0 x 10^^ years is three 
orders of magnitude greater than the predicted 2v half live TqI^ = 1.4 x 10^^ years, which 
essentially agrees with the experimental one, confirming that the observed /3/3 decay oi'^^^U has 
to be the two neutrino mode. Figure 1 exhibits the M^J [jy) and — M(^(ji,) components of the 
Gamow- Teller and Fermi nuclear matrix elements for "^^^U . It can be seen that all the terms 
contribute constructively, and that the transition is strongly dominated by a pair of neutrons in 
the normal parity orbital z^j'j^ ^ decaying into two protons in the unique orbital 'i\-^i2-, resembling 
the two neutrino case (see [jl2| and Appendix B). 

In the case of ^^^Nd, the pseudo SU(3) Ov matrix element reported here is about a factor four 
lesser than the QRPA estimations. This is a very relevant result. First, it exhibits the stability of 
the neutrinoless double beta decay matrix elements evaluated in quite different nuclear models, 
in the case of deformed nuclei. Second, this factor of four, which is little compared with the order 
of magnitude variations in the 2z^ theoretical estimations, is still important in order to extract 
the parameters < rrie > and < g^^y^ >. In Figure 2 the M^J'{jy) and — Mq^(J[,) components 
of the Gamow- Teller and Fermi nuclear matrix elements for ^^^Nd are shown. As for 238^ ^yi 
the terms contribute constructively, but in this case the transition is dominated by a pair of 
neutrons in the normal parity orbital /ig /g decaying into two protons in the unique orbital h'^i/2j 
again resembling the two neutrino case. 

1/2 

As can be seen in the last two columns of Table 1, the Tq^ predicted for < niu >= leV 
are at least three order of magnitude greater than the experimental limits. These results reflect 
the fact that at, the present stage of the experimental /3/3 research, the limits < rrii, >< l.le^ 
obtained by the Heidelberg-Moscow collaboration ||2^ using significative volumes of ultrapure 
^^Ge are the most sensitives. But, if the PPou decay is observed in ^^Ge, at least a second 
observation will be essential, and ^^^Nd is a likely candidate to do this job |16|. In the next few 



years the limit for < m,^ > extracted from PPou experiments is expected to be improved up to 
0.1 eV, and ^^^Nd is one of the selected isotopes |29]. 



On the other side, in extracting the majoron couphng constant the ^^^Nd has set the lowest 

2 

Umit [ll^]. The phase integral for the two emitted electrons and the majoron, with -^ = 1.0, 

dy 

for this nuclei is 

Goi.x = 0-466 X 10~^''yrs~^ (24) 



and the reported lower limit for this decay is [|l6| 

Toil > 5.30 X 10^%rs. (25) 

Using these numbers in Eq. (3) we obtain the limit 

< gu,x >< 4.1 X 10-^ (26) 



This limit is less stringent than the previous one |16| obtained with the same experimental 
information but employing the QRPA matrix elements Mg^ = 6.075 |15]. 

There are various elements which contributes to the greater cancellation of the pseudo SU(3) 
Mojy matrix elements, as compared with the QRPA ones. 

i) The uncertainty in the QRPA results itself: the other QRPA reported value for ^^^Nd is 
Mq^ = 5.24 Q, 14% lesser than the above mentioned one. 

ii) The way short range correlation effects were parametrized in H{r, E) in the present work 
is not the most common one. But both finite nucleon size (FNS) and two-nucleon short range 
correlations (SRC) effects together reduce the Qv matrix elements in less than 15%. 

iii) Another source of cancellation could be related with the seniority zero approximation 
used to describe nucleons in abnormal parity states. Removing this assumption in Eq. ( |T8| ) 
would lead to a more general expression, not restricted to states with J = M = 0. Being far 
from trivial, this extension of the pseudo SU(3) model is under investigation. But as can be 



seen in Fig. 3 of |15], the {j-kJ-k)J = 0; {jvju)J = channel represents by far the most relevant 
contribution to the Qu matrix elements, and inclusion of other multipoles could decrease Mq^. 

iv) Only one active shell was allowed for protons, and one for neutrons. This is a very strong 
truncation, of the same type as used in shell model calculations. Extensions of the pseudo SU(3) 
model which takes fully into account core excitations have shown that the SU(3) component of 
the wave function are the dominant ones, between 50% and 60% of the total amplitude [pQ] . 
Inclusion of these symplectic excitations could increase the present results, allowing for example 
transitions from neutron states to proton states in the same shell. Given the SU(3) dominance, 
these corrections are estimated to contribute at most with a factor two. 

Nevertheless, at the present stage of development, the pseudo SU(3) predictions are as valid 
as the QRPA ones, which suffer other kind of limitations, and their differences are in some 
way describing the theoretical uncertainties in the extraction of the parameters < niu > and 

< 9u,x >• 



6 Conclusions 



Our description of the neutrinoless /3/3 decay in heavy deformed nuclei can be summarized as 
follows: 

The pseudo SU(3) model is a very powerful machinery to describe the collective behavior of 
heavy deformed nuclei. Together with a model hamiltonian with spherical Nilsson single particle 
terms, and a residual interaction of a quadrupole-quadrupole type, plus terms including some 



9 



scalar combinations of the quadrupole and orbital angular momentum operators, it has been 
used to reproduce very accurately the rotational spectra of heavy deformed nuclei, including 
the K-band splitting, the amplitudes for transitions of the E2, Ml and M3 type and the 2v(3(3 
decay of eleven heavy deformed nuclei, with good agreement with the available experimental 
data. 

The deformed ground states were described using the pseudo SU(3) scheme, which is able 
to generate maximum deformation for nucleons in normal parity orbitals, filling them into the 
Nilsson levels but with definite total angular momentum. These states are strongly correlated, 
and are dominated by the proton-neutron quadrupole-quadrupole interaction, considered the 
microscopic origin of deformation. The nucleons in the abnormal parity orbitals are assumed to 
have seniority zero, i.e. in this approximation they do not contribute actively to the quadrupole 
moments. This is the strongest approximation of the model, only justified in terms of its 
simplicity. 

We have evaluated the I3I3qu half lives of six heavy deformed nuclei, using the pseudo SU(3) 
approach together with the closure approximation, obtaining an order of magnitude agreement 
with the QRPA estimations, and a factor four of difference. We exhibited predictions for the 
neutrinoless double beta half-lives assuming < rriy >= leV , and discussed the relevance of our 
results. In the case of "^"^^U this result complements those obtained for the two neutrino I3j3 
decay, and it confirms the observed half-life as two neutrino in origin. Limits for the majoron 
coupling constant were extracted from the ^^^Nd experimental limit, and compared with the 
QRPA ones. 



10 



Appendix A: The radial integrals 



The radial integrals are defined as|25] 



oo 




oo oo 



{/ rfdriji{pri)Rnii^{ri)Rn2i2{r2) J Hdr23i{pr2)Rn-ih{ri)Rn^i^{r2)]dp 



v{p) = f / H{r)jo{pr)r'^dr 


(27) 
The Rni are the one-particle radial harmonic oscillator wave functions and jo{pr) is the 0-th 
order spherical Bessel function. 



Using the explicit form of H(r,E), Eq. (1), we find |31] 



1 

vr p(p -|- E) 



The explicit expression for the radial integral is [^5[ 

R\hl2,hk) = iM{nih,mh)M{n2l2,nAh))-^^^ 

J2 amAriih,n3l3)amAn2l2,n4i)a2m{^^l,^^l)Jm ^'^^^ 

ml,m2,m 

where mi = li -\- 1^ + 2s and s = 0, 1, ...,?ii + n^, and similarly to 1712, m = 1,1 -\- 1, ..., Hli±IM^ 
M{nl, n'l') = 2"+"'n!n'!(2/ + 2n + 1)!!(2/' + 2n' + 1)!! 

n, „ . (r^l r,'l'\ - ( ^Y V n! n'\ {2l+2n+l)\\ {2V +2n' +l)\\ (30) 

fl,fJ,'=S — fl 

and 



00 



Jm = {b^/2r I exp{-^)v{p)p'^+^dp (31) 



being b = \ — the oscillator length. 
In the case £■ = 0, we obtain 

H{r,E) =H{r,0) - ^ 



v{p) = ^ Jm = #r(m + 1/2) 



(32) 

which is a not so bad approximation to the exact result ||17|], obtained using a nonzero E and 
including the short range correlations and finite nucleon size effects, as explained in Appendix 
C. 



Appendix B: The two-neutrino limit 
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It is interesting to exhibit how the two neutrino matrix elements presented in |12] can be 



obtained from the present formahsm. In the two neutrino case H(r) -^ 1, and the transition 
operators are the double Fermi or the double Gamow- Teller ones. In this limit 

/ Jl{pri)jiipr2)v{p)p'^dp= - H{r)Pi{cosuj)d{cosuj) ^ 5io (33) 



where r = \fi — f2\ and oj is the angle between fi and r2 and 

< {lil2)LM\V{r)\{kU)LM > -^ Sn.n-Aih^n^nJhk 

where rij are the number of nodes in the radial wave function Rmk ■ 

Using the previous results and expressing j^^ = 1-k + 1/2, i;/ = I-k — 1/2 we obtain 



(34) 



E(2S + 1)W{^1^ + ^Sl^; l4)Wi^l^ - ^Sl^; Z^) < {^^)S\\T ■ r||(ii)5 > 

= Sn.nJu.^/^4^{< i-2-2)s = o\\T ■ r||(ii)5 = > - < (ii)5 = i||r • r||(ii)5 = i >}. 

(35) 
In the case of double Fermi transitions, i.e. F • F = J2mn ^m^n i the above described matrix 



elements are always zero. This result was argued heuristically in |12], and exhibit the isospin 
as a good quantum number for the pseudo SU(3) ground state. For double Gamow- Teller 
transitions, where F • F = X^mn ^m ■ ^ntm^n ^^ obtain 



< {J.J.)J = OIO^^KjW.) J = >^ -5^^^J,^,^A^^.^^JL±^ (36) 

which substituted in Eq. (18) gives the < Ot|F • F|0,^ > in exact equivalence with Eq. (6.8) 
and (6.9) of |jl2|. This 2v double Gamow- Teller operator creates a pair of protons coupled to 
total angular momentum zero and annihilates two neutrons of the normal parity space coupled 
to pseudo orbital angular momentum and pseudo spin equal to zero. For the rare earth and 
actinide nuclei this operator annihilates two neutrons in the pseudo shells r]^ = 5 and rj^ = 6 
and creates two protons in the abnormal orbits /111/2 and ^13/2, respectively. 



Appendix C: The finite nucleon size and short range correlation effects 



In the momentum representation it is not difficult to include additional effects. The finite 
nucleon size (ENS) effects are introduced through the dipole form factor |^ p^ 

VFNsip) = Hp){ ,2 , 2 }' ^ = 850MeV (37) 

The two-nucleon short range correlations (SRC) mainly arise from the repulsion due to the 
cj-exchange in the nucleon-nucleon interaction ||3^. Thus the correlated potential is constructed 
as 
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n^p.k)^^27TfS{p)-^5{q-qc) ^ ^ 

with Qc = 3.93/?7i~^ roughly the Compton wavelength of the a; meson. Within the approxima- 
tion Qc ^ E we get 

vsRcip) = v{p) - Av{p) Av{p) = — ln\ _ "" I (39) 

QQc Q Qc 

When both the FNS and the SRC effects are considered the neutrino potential takes the 
form: 

VFNS+SRC{P) = VFNS " Af (p) + Av' (p) (40) 

with 



M _ . , . x" . „ A2 



Av'ip) = ^y -[x" - xl] + In—} xl = — -^. (41) 

Pqc^,n^~ "^ xl^ ^ A2 + (p±g,)2 ^ ^ 



It was shown in Ref. |31| that the FNS effects are negligible for the allowed transitions but 
increases rapidly when L increases, and that, if both FNS and SRC are considered simultaneously 
the net result may not be very different from that obtained without any SRC at all. 
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Table Captions 



Table 1. Theoretical estimates for the nuclear matrix elements and the half-life of the /3/3-decay 
in the Ov mode for several heavy deformed nuclei are given, under the assumption < me >= leV, 
and compared with the available experimental data. 



Table 1 



Transition 


Mou 


Gouiyrs 1] 


E[MeV] 


^1/2 
^theo 


Texp 


146 iVd - 


-^146 Sm 


1.16 


1.65 X 10-17 


13.5 


1.18 X 1028 




148 iVd - 


-^148 Sm 


1.57 


1.57 X 10-14 


13.6 


6.75 X 1024 




150^^ - 


-^150 Sm 


1.57 


1.01 X 10-13 


13.7 


1.05 X 1024 


> 2.1 X 102i ll^ 


186 p^ _ 


_^186 Q^ 


1.70 


1.76 X 10-15 


15.3 


5.13 X 1025 


> 2.3 X 1020 |33| 


1920s - 


-^192 Pi 


0.72 


1.58 X 10-15 


15.5 


3.28 X 1026 




238 f/_ 


-.238 pu 


1.76 


8.21 X 10-14 


17.3 


1.03 X 1024 


> 2.0 X 1021 
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Figure Captions 



Figure 1. The MqJ'{Ji,) (white columns) and — Mq^(Ji^) (black columns) components of the 
Gamow- Teller and Fermi nuclear matrix elements for ^^^[/, as function of the angular momentum 
of the decaying neutrons. 



Figure 2. The same of Fig. 1 for '^^"Nd. 
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Figure 1: 



Figure 2: 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9407022vl 



This figure "figl-2.png" is available in "png" format from: 
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